Calculated Magnetic Exchange Interactions in High- Temperature Superconductors 
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Using a first principles linear response approach, we study the magnetic exchange interactions J 
for a series of superconducting cuprates. We reproduce the observed spin-wave dispersions together 
with other experimental trends, and show that different cuprates have similar J's regardless their Tc. 
The nearest neighbor J is not sensitive to the hole-doping, which agrees with recent experiments. 
For the undoped cuprates, the second nearest neighbor J is ferromagnetic, but changes its sign with 
hole-doping. We also find that, in contrast to the hopping integral, the exchange interaction is not 
sensitive to the position of apical oxygen. To see the effect of the long-range nature of the exchange 
on the superconducting Tc, we study the dynamical spin susceptibility within the t-J model 

using a dynamical cluster approximation. 

PACS numbers: 75.30.Et, 71.15.-m, 75.10.-b, 74.20. Mn 
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INTRODUCTION 



After more than two decades of intensive studies of 
high-temperature superconductors (HTSCs), it is now 
well accepted that magnetic exchange interactions in the 
cuprates play a fundamental role. Despite a general lack 
of consensus on the pairing mechanism in the cuprates, 
it is believed to be of magnetic origin fl], 0, 0, 0,1^1 • In 
support of this scenario, a direct relationship Tc (x J be- 
tween Tc and the in-plane exchange interaction J was 
extracted from magnetic measurements for one family 
of the cuprates in Ref. |6j|. Despite vast efforts devoted 
to understanding the magnetic properties of the cuprates 
there are still several important issues which need 
to be clarified. Among them is the long range nature of 
the exchange interactions J between the Cu-ions, and its 
effect on the dynamical spin susceptibility x('Z, i-^) and the 
critical temperature Tc. Studies of model Hamiltonians, 
such as the t-J model show a strong dependence of the 
calculated properties on the value of J 22, 2^. In partic- 
ular, numerical calculations suggest a direct relationship 
between the magnitude of the r nag netic exchange cou- 
pling and the pair binding energy i23j. Therefore, obtain- 
ing accurate estimates for the exchange constants and 
studying their trends across different HTSCs including 
their doping dependence and influence on ^(q, lo) and Tc 
is an interesting problem which we address in the present 
work. 



Experiments provide estimates for the magnitude of 
the exchange interaction, but despite many studies there 
is a spread in obtained values even for the nearest- 
neighbor magnetic coupling Ji in the most studied com- 
ound La2Cu04 where it varies from 110 to 150 meV 
ll| . The latter result does not only depend on the 



der exchange terms are included. Nevertheless, based on 
recent high-resolution inelastic neutron scattering (INS) 
experiments, there is an overall agreement that at least 
for the undoped La2Cu04[lH, a simple nearest-neighbor 
Heisenberg model is not sufficient and the high-order 
magnetic exchange interaction is quite important. How- 
ever, all high-order terms have similar effects on the 
spin-wave dispersion and its intensity dependence, there- 
fore, even INS, the most powerful technique for exploring 
magnetic excitations, cannot determine the magnetic ex- 
change coupling very accurately. Worse than La2Cu04, 
for most other HTSC compounds, the difficulty in syn- 
thesizing large single crystals limits the capability of per- 
forming accurate measurements. Thus, the experimental 
information on the exchange interactions is very limited. 

In addition to the experiment, the magnetic interac- 
tions of the undoped cuprates have been studied in a 
number of theoretical works figj. 25|,[2^, 27|. However, 



experimental technique being used but also on how the 
results are fitted, i.e., depending on whether higher or- 



most of them use a cluster approximation and map the 
total energy differences obtained by a first principle cal- 
culation to a model Hamiltonian which prevents detailed 
studies of the long range nature of the J's. It is, for ex- 
ample, well known that the long-range hopping integrals 
are very important for the electronic properties in the 
cuprates [2^ and even the effective third nearest neighbor 
hopping t" has a considerable effect on the interatomic 
exchange [13. Therefore it is desirable to re- investigate 
the exchange interaction in cuprates by a first principle 
calculation. 

In contrast to the undoped cuprates, investigations of 
J's in doped HTSCs are scarce, regardless the supercon- 
ductivity actually happens only after introducing doping. 
Based on a model calculation. Si et at, suggested that 
Ji will decrease rapidly with doping [3l|. On the other 
hand, the commonly used t-J model uses the same J for 
different doping levels [2^. therefore study the effect of 
doping on the exchange interaction based on a first prin- 



ciple calculation is quite interesting. Unfortunately, ab 
initio techniques for extracting J's (like the cluster ap- 
proximations mentioned above) are not efficient for the 
doped case, and there are no extensive calculations of 
this type reported in the literature. In this work, we 
use a recently developed linear response approach [s^ . 
and perform detailed studies of exchange interactions for 
both parent and doped HTSCs. Beyond this, we also 
investigate the effect of the long-range J's on the spin 
susceptibility xil, and the superconducting transition 
temperature Tc within the framework of the dynamical 
cluster approximation (DCA) for the t-J model. 

EXCHANGE INTERACTION AND SPIN WAVES 
Method 



2 




FIG. 1: Definitions of nearest neighbor, Ji, next nearest 
neighbor, J2, and third nearest neighbor, J3, exchange inter- 
actions for Cu spins, which are the parameters of a Heisenberg 
model H = Y.{ij) JijSiSj. 



We perform our electronic structure calculations based 
on density functional theory (DFT) within the full po- 
tential linearized-muffin-tin-orbital (LMTO) method [33l|. 
To take into account the effect of the on-site electron- 
electron interaction we supplement the local density ap- 
proximation (LDA) to DFT by adding a correction due 
to Hubbard U using so-called LDA-I-U approach ^S^], 
with the parameters U=10 eV and J=1.20 eV for the 
Cu d orbitals as deduced from the constrained LDA 
calculation (3^. Experimental lattice parameters have 
been used for all materials. 

With the electronic structure information, one can 
evaluate the magnetic exchange parameters J of a 
Heisenberg model H — J^ij JijSi • Sj based on a mag- 
netic force theoremjs^ which assumes a rigid rotation of 
atomic spin. In this formalism, the interatomic exchange 
constant J is given as a second derivative of the total 
energy difference induced by the rotation of moments at 
sites R + T and R + t [32] : 

X {i^^^^'\[^xB^,]p\^^,)e^^<''~'^\ (1) 

where /, cr, and B are the Fermi function, Pauli ma- 
trix, and the effective magnetic field in the calculation 
given by the difference in the electronic self-energies for 
spin-up and spin-down electrons, respectively. "0 and e 
are the eigenstate and eigenvalue from the LDA-I-U cal- 
culation. This technique has been used successfully for 
evaluating magnetic interactions in a series of transition 
metal oxides[32|- 

Parent Compounds 

We now discuss our results for La2Cu04 and find that 
the exchange constants Ji, J2 and J3 within the Cu02 



plane (see Fig. 1) decrease rapidly with increasing dis- 
tance between two Cu ions. After taking into account 
the effect of quantum renormalization 24], the measure- 
ment of two-magnon Raman scattering gives Ji = 116 
meVQ, which is slightly larger than our numerical Ji 
(109 meV) as shown in Table I. Early neutron scattering 
experiments give a larger value of Ji ^^130 meV0[ld|, 
which may be partially due to the use of only a nearest- 
neighbor Heisenberg model to fit their spin-wave veloc- 
ities. Our Ji agrees very well with one cluster calcu- 
lation (105 meV)[25] but is smaller than the result of 
other calculations (^140 mcV) 19, 26]. Turning to the 
discussion of the next nearest neighbor coupling, our cal- 
culated J2 is slightly larger than the one deduced from 
neutron scattering[9]. In contrast to the early cluster 
calculations fi^j. our J2 is ferromagnetic (FM) and thus 
enhances the antiferromagnetic correlations. 

Since well-defined spin-wave excitations throughout 
the Brillouin Zone have been observed by the INS[ll|, 
it is interesting to perform the comparisons with our cal- 
culated spin-wave dispersions. For cuprates the quantum 
fluctuations may be large due to the smallness of the spin 
S=^ and the low dimensionality D=2. On the basis of the 
Holstein-Primakoff transformation using 1/S expansion, 
it has been found that a renormalization factor is neces- 
sary for the spin-wave excitation energy ^23] in order to 
compare it with the result of standard linear spin-wave 
theory [23]. With this correction, the spin- wave disper- 
sion can be expressed as: 

E, = 2Z,.^A\ - Bl 

where Aq = J\ — Ji\V — cos(27rga;) cos(27rgy)] — J-i\V — 
i(cos(47rgj;) -|- cos(47r(7y))], = i Ji [cos(27rg:r) + 

cos(27rgy)], and Zc is the renormalization factor, respec- 
tively. Here based on the above formula with the quan- 
tum renormalization factor Zc=1.18[2j], and using the 
obtained numerical Ji , J2 , and J3 , we calculate the spin- 




(3/4,1/4) (1/2,1/2) (1/2,0) (3/4,1/4) (1,0) (1/2,0) 

FIG. 2: Comparison between calculated (solid lines) and 
experimental (symbols) spin- wave dispersions for La2Cu04. 
The triangles and squares are the experimential results at 
T=10 K and 295 K, respectively [ij]. 
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FIG. 3: Calculated spin-wave dispersions for 

HgBa2CaCu206, Sr2Cu02Cl2 and YBaaCusOe. 



wave dispersion for La2Cu04, and display the result in 
Fig. 2 by the sohd line. For coinparison, in Fig. 2 we also 
show the INS data by symbolsfllll. Our results agree well 
with the experiments near the zone center. A small dis- 
crepancy exists around the zone boundary which may be 
due to the four-particle cyclic exchange interaction J^., 
which recently attracted much attentionflTl, . 

We now turn to our predictions of exchange interac- 
tions for other HTSC materials where the experimental 
values of J are limited. The experiment [14] |.15j shows 
the Ji of Sr2Cu02Cl2 is about 10 meV smaller than that 
of La2Cu04. This is reproduced by our numerical cal- 
culation as shown in Table I. The Ji of YBa2Cu3 06.i5 
has been measured by neutron scattering, and the ob- 
tained value is about 120 meV[ll| 1131 1. After considering 
the quantum renormalization effect j24|, the experimen- 
tal Ji will be reduced to 100 meV, which is very close 
to our numerical result 93 meV. It is interesting to note 
that our theoretical study reproduces the experimental 
trend across the materials studied (La2Cu04 has largest 
Ji, Sr2Cu02Cl2 is intermediate, and Ji in YBa2Cu306 is 
smallest (lif.) As it is seen from Table I, all parent HTSCs 
have almost the same Ji (around 110 meV) regardless 
of their number of Cu02 layers and their different T^'s. 
J2 is also similar and shows FM behavior while J 3 is 
AFM like. Using the quantum renormalization factor 24 
and linear spin-wave theory, we also calculate the spin- 
wave dispersion for all other compounds and show the re- 
sults for HgBa2CaCu206, Sr2Cu02Cl2 and YBa2Cu306 
in Fig. 3. Since different compound have similar exchange 
interactions the shape of the spin-wave curve is quite sim- 
ilar while YBa2Cu306 has smaller spin-wave excitation. 

To further check the possible relationship be- 
tween J and Tc, we also study a high-pressure 



phase of H gB a2 Ca2Cu30s, for which Tc increases 
considerably [37|. Applying pressure reduces the lattice 
parameter, and, as a result, increases the magnitude of 
J, but this enhancement is too weak to relate it with ris- 
ing Tc. For example, we find the values of Ji are 114.3 
and 116.4 meV at 3 and 8 GPa, respectively, which is 
very close to the value at ambient pressure as shown in 
Table 1. 



The effect of apical oxygen 



It is believed that the apical oxygen has a dramatic 
effect on TcjH, 4^ 41, i^]. However, there is a debate on 
whether Tc is positively correlated with the ma gnit ude 
of the effective next nearest neighbor hopping t' [4111 , or 
if an intersite super- repulsion term is important |42[ | . We 
therefore study the effect of the position of the apical 
oxygen on the exchange interaction. The apical oxygen in 
La2Cu04 is located at the (0,0, Zq) site. By adjusting the 
internal atomic coordinate Zq we perform the calculation 
for different distances d-A between the apical oxygen and 
Cu. 



Our results are shown in Table II. In contrast to the 
hopping integral, which is quite sensitive to the position 



of the apical oxygen [4 iL [42[ |. the magnitude of the ex- 



change interaction in the Cu02 plane only slightly in- 
creases with d-A as shown in Table II. In the Hg-based 
cuprates with d^ ~ 2.7 A J is found smaller than the 
one in La2Cu04 with the same d^, so it is likely that the 
exchange interaction is more sensitive to the detail of the 
electronic structure, rather than just to the distance d^. 
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TABLE L Experimental Tc (K), and calculated exchange in- 
teractions (meV) for parent HTSC materials. Ji, J2, and 
J3 are the nearest-neighbor, second nearest-neighbor and 
the third nearest neighbor exchange interactions as shown in 
Fig.l. Njai/er is the number of CUO2 layer, Tc is critical tem- 
perature. 



Jl 



J2 J3 



CaCu02 


1 




110.0 


-10.1 


3.8 


Tl2Ba2Cu06 


1 


97 


109.1 


-10.9 


4.0 


HgBa2Cu04 


1 


94 


108.9 


-11.1 


3.3 


La2Cu04 


1 


42 


108.8 


-12.0 


-0.2 


Sr2Cu02Cl2 


1 


28 


99.2 


-8.2 


1.6 


HgBa2CaCu2 06 


2 


128 


110.4 


-11.9 


2.9 


Tl2Ba2Cu208 


2 


125 


108.7 


-10.7 


2.5 


YBa2Cu306 


2 


90 


93.0 


-4.7 


2.4 


HgBa2Ca2Cu308 


3 


135 


109.9 


-10.1 


2.8 



TABLE II: The calculated exchange interaction in La2Cu04, 
with different A a, where A a is the distance between apical 
oxygen and Cu atom. Aa is in A and J is in meV. 



dA 


Jl 


J2 


J3 


2.5 


111.1 


-12.4 


-0.4 


2.6 


112.9 


-13.1 


0.1 


2.7 


114.2 


-13.8 


1.2 


2.8 


116.0 


-14.6 


2.1 



The effect of doping 



It may not be surprising that the J's of the par- 
ent HTSCs do not directly correlate with their Tc's 
which characterize the corresponding doped materials. 
We therefore study the effect of hole doping for La2Cu04 
using the virtual-crystal approximation (VGA), which 
has been used successfully for the phonon properties of 
La2-2:Ba:rCu04[3l]. For the doped case, our scheme is 
rough, but nevertheless does include major ingredients 
of the system, such as, superexchange, double exchange 
and RKKY exchange interactions 3l| . Naively, one may 
think the hole induced by doping can hop through the Cu 
ion, and result in a ferromagnetic like double-exchange 
interaction, which would consequently suppress Ji. But 
if the doping level is not high, this effect is not large as 
seen in Table III. Since the spin-wave velocity is mainly 
controlled by Ji , our results agree with the recent INS ex- 
periments, which showed clearly that for La2-a;Sra;Cu04 
the spin-wave velocity is doping insensitive [l6!|. Hole 
doping enhances J3 slightly as seen in Table III. Differ- 
ent from Jl and J3, doping has a large affect on J2, which 
changes from FM-like to AFM-like resulting in consid- 
erable spin fluctuation. 



TABLE III: Doping effect on exchange interactions in 
La2Cu04, where x is the hole-doping concentration. J is 
in meV. 



x 


Jl 


J2 


J3 


0.0 


108.8 


-12.0 


-0.2 


0.1 


110.9 


-7.9 


-3.0 


0.2 


117.8 


-0.3 


-3.1 


0.3 


124.6 


6.4 


-3.8 



EFFECT ON Tc 



Now, we aim to address the following questions: (1) 
Can the variation in exchange parameters for different 
materials explain the difference in Tc? (2) What is the 
effect of varying exchange parameters on the dynamic 
spin fluctuation spectrum? To answer these questions, 
we use a dynamic cluster approximation (DCA) [il, Q 
to calculate the properties of a t- J model. A similar study 
was performed in Ref. [i^. There, a combined DFT- 
LDA and DCA-QMC approach was used to study the 
parameter dependence of Tc in a three-band Hubbard 
model. These calculations showed that Tc is a very strong 
function of the hopping parameters, and very sensitive to 
even small variations in the long-range hopping integrals. 

Here, we want to study the dependence of Tc on the 
exchange parameters J. We therefore consider a two- 
dimensional t-J model 



where Si = cj cr. 



+ c ■ ) 

IS J s I 



E 



Jij Si Sj , (2) 



and cj„ 



is'-' ss' ^is' "ii^ ^is is projected fermion op- 
erator deflned as c^s(l — rii^s)- 

The general idea of the DCA is to map the bulk lat- 
tice problem onto an effective periodic cluster embedded 
in a self-consistent dynamic host designed to represent 
the remaining degrees of freedom. Correlations within 
the cluster are treated explicitly, while those beyond the 
cluster size are treated on the mean-field level. The hy- 
bridization of the cluster to the host accounts for fluctu- 
ations arising from the coupling between the cluster and 
the rest of the system. 

The mean-field nature of the approach allows us to 
study transitions to symmetry broken phases such as the 
superconducting state even in small clusters. For exam- 
ple, using a dynamic cluster quantum Monte Carlo ap- 
proximation for a small four-site 2x2 cluster, the proper- 
ties of a 2D Hubbard model were calculated in Ref. [46|. 
The obtained phase-diagram is remarkably similar to the 
universal cuprate phase diagram, exhibiting antiferro- 
magnetic and d-wave superconducting phases as well as 
pseudogap behavior. As discussed in Ref. ^il], a four-site 
cluster DCA calculation provides a mean-field result for 
the transition to a superconducting state with a dx2_y2- 
wave order parameter. DCA calculations for larger clus- 
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FIG. 4: Spin susceptibility x'(q,uj) for q = (7r,7r) versus u) for 
different exchange interactions J2 and J3. Here, Ji = 0.3t is 
fixed and a filling (n) = 0.8 and temperature T = O.OSf have 
been used. 
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FIG. 5: Superconducting transition temperature for different 
exchange interactions J2 and J3. Here, we have used the same 
values for Ji , J2 and J3 and the filling as in Fig. |4] 
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ter sizes show a reduction in the superconducting due 
to the inclusion of pair- field phase fluctuations, but the 
qualitative aspects of the results including the pairing 
mechanism are similar to the four-site cluster results [4^ . 

Here, we have used a non-crossing approximation 
IS, H [S^l (NCA) to determine the spin susceptibility 
x(g, w) and Tc for a four-site 2x2 cluster. We perform 
the simulations in the superconducting state by allow- 
ing for a finite anomalous Green's function With 
increasing temperature, Tc is determined by the temper- 
ature where the the anomalous Green's function vanishes. 
Similar calculations for fixed near-neighbor exchange in- 
tegral J = 0.3i were performed in Ref. [H^l- Many results 
of this study were shown to be reminiscent of experi- 
ments in the cuprates- Here, we focus on the effect of 
varying longer-ranged exchange parameters on the spin- 
susceptibility and the superconducting Tc- 

Results showing the imaginary part of the spin suscep- 
tibility, x"((7,a;) versus uj for q — (7r,7r), calculated at a 
temperature T — 0-08i and fixed Ji = 0-3t, are plotted 
in Fig. [5] for various values of the next-nearest-neighbor 
and third-nearest neighbor exchange interactions J2 and 
J3. As one can see, a ferromagnetic J2 < and an an- 
tiferromagnetic J3 > enhance the spectral weight in 
x"{q,Lo) at the antiferromagnetic wave vector q = (7r,7r) 
at low frequencies. 

Previous DCA quantum Monte Carlo and NCA simu- 
lations have addressed the question of the pairing mech- 
anism in the Hubbard and t-J models [13]- The results 
of a recent DCA/NCA and Lanczos study of the super- 
conducting gap function were found to be consistent with 
a simple phenomenological form for the c?-wave pairing 



Vd{k,L0,k\L0') = (3) 

3 - 

-f7^x(fc — k',uj — u') — J(cosfc2; — cosA:j^)(cosA;^ — cosk'y) . 

Here, tJ and J are effective coupling constants of the 
retarded spin fluctuation part and the non-retarded ex- 
change contribution, respectively. It was shown that 
the dominant contribution to the pairing interaction Vd 
comes from the spin fluctuations [Sli] . 

From this one would expect that variations in the spec- 
tral weight in ^"(g, cj) will directly affect the strength of 
Vd and thus T^. Fig. [5] shows the DCA/NCA results 
for Tc for different values of the next-nearest-neighbor 
and third-nearest-neighbor exchange integrals J2 and J3. 
Here, we have fixed the nearest-neighbor exchange Ji = 
0.3i and the filling (n) = 0.8. Consistent with Eq. dS]), 
one finds that the magnitude of Tc tracks the magnitude 
of the spectral weight in x" {1^ ^) for q = (tt, it). E.g., the 
highest Tc is obtained for J2 = — 0.2Ji, J3 — 0, i.e. when 
the spectral weight in x"((7,a;) is maximal for q = {tt,tt). 

We emphasize, however, that is rather insensitive to 
changes in the long-range exchange parameters J2 and 
J3. E.g. a change in J2 from J2 = to J2 = — 0.2Ji only 
induces a ~ 5% increase in Tc and the effects of J3 on 
Tc are almost negligible. Note, however, that the third- 
nearest neighbor exchange J3 is not directly represented 
on a 4-site 2x2 cluster. In this approximation, the spins 
on the cluster only know about J3 through their effec- 
tive coupling to the spin-fluctuation host. Therefore, it 
is possible that larger variations in Tc would be found 
in larger cluster simulations. ^From the present results, 
however, we conclude that the differences in the long- 
ranged exchange terms cannot explain the differences in 
Tc between different cuprates. 
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CONCLUSION 

In summary, we have calculated the magnetic exchange 
interactions for a series of HTSC cupratcs using the 
LDA+U based linear response approach. Our calculated 
spin wave dispersions were found in good agreement with 
high-resolution inelastic neutron scattering experiments. 
Our simulations show that for different parent HTSCs, 
the nearest neighbor and next nearest neighbor exchange 
constants are similar, while the third nearest neighbor J 
appears to be material dependent. Furthermore, we find 
that Ji is insensitive to the hole doping, which agrees 
with recent experiments and supports the implicit as- 
sumption of the t-J model with a fixed magnitude of 
J. In contrast to the hopping integral, which strongly 
depends on the position of the apical oxygen, our results 
show that the magnetic exchange interaction is rather in- 
sensitive to the position of the apical oxygen. In addtion 
to the LDA-I-U calculations, we studied the dependence 
of the dynamic spin susceptibility x(9, and supercon- 
ducting transition temperature Tc on the exchange pa- 
rameters in a t-J model using a dynamic cluster approxi- 
mation. We find that both x{q, and Tc are only weakly 
affected by variations in the exchange parameters beyond 
nearest neighbor. Based on these results, wc; conclude 
that differences in long-range exchange terms between 
different materials cannot explain their different super- 
conducting transition temperatures. 
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